<1*727667 


RR  288 


Research  Report  288 


INFLUENCE  OF  SOIL  ON  DETECTION  OF 
BURIED  EXPLOSIVES  AND  TUNNELS 
BY  TRACE  GAS  ANALYSIS 

R.P.  Murrmann.Y.  Nakono, 

T.J.  Simpson,  D.C.  Ltgoott  and  D.M.  Andorson 


Juno  1971 


Roproducod  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Spnngftold,  Va. 


CORPS  OF  ENGINEERS,  U.S.  ARMY 


DDCr 

pnn.njE 

MG  * 

HSEUTj'E 


COLD  REGIONS  RESEARCH  AND  ENGINEERING  LABORATORY 


HANOVER,  NEW  HAMPSHIRE 


43 


APPROVED  FOR  PUBLIC  RELEASE!  DISTRIBUTION  UNLIMITED. 


Unclassified 


37 


REF’ORt  SECURITY  C  L  ASSIMC  A  tiQK 


Unc lassi f ied 


Sfctmlv  Has^ific  ifion 


DOCUMENT  CONTROL  DATA  R  &  D 

v**i  •  '*  ft  ■  >h>'  i/i-  nfion  of  title.  Inulj  of  mhitrt n  t  anil  indexing  annotation  must  be  entered  nfirn  ffii*  overall  report  l  %  < 


U^S.  *Xrmy*  dolid’  fee^ions^ftesearch  and 
Engineering  Laboratory 
Hanover,  New  Hampshire  03755 


Kl  Cl'  I  t|  L  k 

INFLUENCE  OF  SOIL  ON  DETECTION  OF  BURIED  EXPLOSIVES  AND  TUNNELS 
BY  TRACE  GAS  ANALYSIS 


4  PEiCRiPTivE  notes  (T'-pa  ol  report  and  Inclusive  date*) 


■>  au  t  ho  RiS»  fA-'irif  n#m«.  middle  Initial,  laal  name) 

R.P.  Murrmann;  Y.  Nakano;  T.J.  Simpson;  D.C.  Leggett  and  D.M.  Anderson 


6  REPOR T  DATE 

June  1971 


l«.  CONTRACT  OR  GRANT  NO 


6.  PROJIC  T  NO 


7R.  TOTAL  NO  OF  PACKS  1 7b.  NO  OF  REF* 


9a.  ORIGINATOR’S  REPORT  NUMBER!*) 

Research  Report  288 


DA  Project  1J662708A462 


9b.  OTHER  REPORT  NOIII  (4ny  other  numbmra  that  aiay  be  aaalgnai 
thla  rapott) 


10  DISTRIBUTION  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


12.  SPONSORING  MILI  T  ARY  ACTIVITY 


U.S.  Army  Mobility  Equipment  Research 
and  Development  Center 
Ft.  Bel voir,  Vu. 


I*  ABSTRACT 


Detection  of  mines,  explosives,  and  tunnels  may  be  accomplished  by  sensing 
associated  volatile  effluvia.  This  investigation  was  undertaken  to  provide 
a  basis  for  predicting  the  diffusion  of  volatile  compounds  from  underground 
sources  into  the  atmosphere.  Diffusion  of  a  volatile  compound  was  studied 
for  a  range  of  soil  conditions  utilizing  soils  from  the  mine  detection  sites 
in  Puerto  Rico.  A  new  mathematical  analysis  based  on  the  Monte  Carlo  method 
was  developed  for  predicting  vapor  diffusion  through  soil  into  the  atmosphere. 
It  was  determined  that  diffusion  in  soil  can  be  reliably  predicted  if  soil 
porosity,  moisture  content,  and  affinity  for  the  compound  are  known.  Appear¬ 
ance  in  the  atmosphere  is  also  dependent  on  accumulation  of  the  compound  in 
air  at  the  soil/atmosphere  interface.  Diffusion  of  volatile  compounds  through 
soil  into  the  atmosphere  is  not  likely  to  be  an  important  factor  in  tunnel 
detection  due  to  depth  of  overburden.  However,  adsorption  of  compounds  at 
tunnel  walls  is  likely  to  significantly  reduce  the  amount  of  vapor  appearing 
in  the  atmosphere  through  entrances  and  vents.  Detection  in  the  atmosphere  of 
TNT  vapor  from  mines  and  explosives  buried  in  moist,  porous  soil  should  be 
possible  under  ideal  sample  collection  conditions. 
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INFLUENCE  OF  SOIL  ON  DETECTION  OF  BURIED  EXPLOSIVES  AND 
TUNNELS  BY  TRACE  GAS  ANALYSIS 

by 

R.P.  Murrmann;  Y.  Nakano;  T.J.  Simpson;  D.C.  Leggett  and  D.M.  Anderson 


INTRODUCTION 

Detection  of  volatile  effluvia  has  been  proposed*  as  the  basis  of  a  method  for  locating 
personnel.  The  same  principle  can  be  applied  to  detecting  underground  discontinuities  such  us 
explosives  or  mines.1  The  applicability  of  this  approach  to  tunnel  detection  has  been  demonstrated.16 
Although  the  chemicals  or  combinations  of  chemicals  upon  which  to  base  such  detection  schemes 
are  still  undetermined,  vapors  originating  from  explosives  or  explosive  containers  could  provide  the 
key  to  detection  of  mines. 

Volatile  chemicals  unique  to  a  tunnel  may  originate  as  a  result  of  perturbation  of  biological 
activity  near  tunnel  walls;  chemical  effluvia  associated  with  human  activities  may  appear  in  the 
atmosphere  over  occupied  tunnels.  In  any  case,  volatile  compounds  originating  from  underground 
sources  can  migrate  into  the  atmosphere  by  diffusion  upward  through  the  soil  profile.  In  tunnels 
they  may  also  move  through  entrances  and  vents.  However,  the  presence  of  soil  between  source 
and  detection  device  complicates  detection.  This  is  because  chemicals  originating  underground 
interact  with  the  soil  constituents,  thus  delaying  their  appearance  in  the  atmosphere  and  decreasing 
the  amounts  available  for  detection  at  the  surface.  Also,  the  relative  amounts  of  chemical  com¬ 
ponents  which  appear  in  the  atmosphere  may  differ  from  those  at  the  undergiound  source. 

This  report  is  concerned  with  predicting  the  extent  to  which  the  presence  of  soil  will  interfere 
with  detection  of  buried  explosives  and  tunnels  by  trace  gas  analysis.  The  effect  of  soil  texture, 
porosity,  and  moisture  content  on  vapor  diffusion  in  soils  collected  from  the  mine/tunnel  detection 
research  sites  in  Puerto  Rico21  22  is  evaluated  in  the  laboratory.  The  extent  of  interaction  of  vapor 
with  mineral  surfaces  in  soil  collected  from  the  Laguna  Joyuda  tunnel  site  is  determined.  A 
mathematical  model  to  describe  diffusion  of  volatile  chemicals  through  earth  media  is  formulated 
which  employs  assumptions  believed  to  be  reasonable  for  application  of  the  model  to  the  problem  of 
detection.  Another  theoretical  approach  based  upon  the  Monte  Carlo  method  is  developed  for  analysis 
of  diffusion  under  different  environmental  conditions.  Calculations  are  made  of  the  time/concentra¬ 
tion  relationships  for  diffusion  of  vapor  through  soil  at  the  Laguna  Joyuda  site.  Parameters  are 
ider-ified  which  are  most  important  in  influencing  diffusion  of  a  volatile  component  through  soil 
into  the  atmosphere  so  that  favorable  environmental  conditions  for  detection  by  trace  gas  analysis 
can  be  anticipated.  Detection  in  the  atmosphere  of  a  volatile  chemical  placed  in  the  soil  at  the 
Laguna  Joyuda  site  is  accomplished. 


♦References  3,  4,  ft,  7.  12,  13,  14,  19,  25. 
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ADSORPTION  OF  VOLATILE  CHEMICALS  BY  SOIL 


Introduction 

The  movement  of  volatile  chemicals  from  subsurface  sources  through  soil  into  the  atmosphere 
is  influenced  by  the  interaction  of  a  given  chemical  with  soil  minerals  and  water.  In  order  to  assess 
the  degree  of  interaction  of  volatile  chemicals  with  soil,  a  simple  method  for  determining  adsorption 
isotherms  lias  been  developed.  The  adsorption  of  acetone,  a  volatile  human  effluvium,  was  deter¬ 
mined  on  soil  collected  from  the  Laguna  Joyuda  detection  site  in  Puerto  Rico.  The  data  obtained 
were  later  used  in  ot her  experimental  and  theoretical  studies  on  predicting  the  effect  of  soil  on 
detection  by  trace  gas  analysis. 


Experimental  procedure 

Soil  for  tliis  experiment  was  collected  from  the  Laguna  Joyuda  site  in  Puerto  Rico.  A 
description  of  the  site  and  the  chemical  properties  of  the  soil  has  been  published  elsewhere.11  11 
Some  of  the  soil  properties  are  given  in  Table  I.  Samples  of  the  soil  were  prepared  with  field- 
moisture,  air-dry  and  oven-dry  water  contents  corresponding  to  about  0.2,  0.02  and  less  than 
0.001  gram  of  water  per  gram  of  dry  soil,  respectively.  About  8  g  of  soil  at  field  moisture  was 
weighed  into  each  of  six  500-ml  Erlenmeyer  Uasks  fitted  with  rubber  stoppers  containing  short  pieces 
of  glass  tubing  flared  at  one  end  to  accommodate  a  silicone  rubber  septum.  The  actual  volumes  of 
the  flasks  ranged  from  545  cm1  to  555  cm!.  Acetone  vapor,  obtained  from  a  similarly  fitted  flask 
containing  liquid  acetone  in  equilibrium  with  acetone  vapor  at  25°C,  was  introduced  through  the 
septum  into  each  flask  by  means  of  a  gas  syringe. 

Table  I.  Description  of  soils  collected  from  detection  sites  in  Puerto  Rico 
for  investigation  of  acetone  diffusion.21  21 


Site  name 

Soil  type 

Mineralogy 

pH 

CEC 

Rio  On  ay  an  es 

Yivi  sandy  loam 

5.0 

17 

S.ibana  lloyos 

Bayainon  sandy  clay 

Kaolinite 

4.2 

14 

Lacuna  Joyuda 

Nipe  clay 

Kaolinite 

5.0 

17 

Lajas 

Aguirre  clay 

Kaolinite  and 
montmorillonite 

6.3 

50 

Tom  Negro 

Cialitos  clay 

Kaolinite 

4.3 

35 

Moca 

Coloso  clay 

Kaolinite  and 
montmorillonite 

6.0 

50 

Measured  volumes  of  acetone  vapor  ranging  from  100  to  1600  pi  were  injected  into  the  sample  flasks. 
A  reference  bottle  was  prepared  by  injecting  1000  gl  of  acetone  vapor  into  a  flask  containing  no 
soil.  An  equilibration  period  of  24  hours  at  25°C  was  allowed  to  pass  before  the  vapor  over  each 
sample  was  analyzed  for  acetone  by  gas  chromatography.  Preliminary  work  showed  that  95%  of  the 
acetone  eventually  adsorbed  by  the  soil  was  taken  up  during  the  first  2.5  hours. 

Following  injection  into  the  chromatograph  of  100  gl  of  vapor  from  each  of  the  flasks,  including 
the  reference  flask,  the  area  under  the  acetone  peak  was  determined  using  a  disk  integrator.  The 
amount  of  acetone  present  in  each  sample  was  determined  by  referring  the  integrated  areas  for  each 
sample  to  a  curve  prepared  from  acetone  in  water  standards.  This  resulted  in  values  for  the  quantity 
of  acetone  vapor  in  equilibrium  with  the  soil  samples.  The  amount  of  acetone  adsorbed  by  the  soil 
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was  then  calculated  as  the  difference  between  the  total  amount  of  acetone  added  to  each  tlask  and 
that  remaining  in  the  vapor  phase.  The  same  procedure  was  followed  using  samples  at  air-dry  and 
oven-dry  moisture  contents  except  that  smaller  quantities  of  oven-dry  soil  (0.5  g)  were  used  per 
sample  because  of  greatly  enhanced  acetone  adsorption  upon  removal  of  the  soil  water.  The 
operating  conditions  for  the  chromatographic  analysis  were  as  follows: 

Instrument:  Perkin-Elmer  Model  900 

Column:  6  ft  *  in.  OD,  SS,  5%  polyethylene  glycol  (1300-1600)  on  40-60 

mesh  chromosorb  T 
Carrier  gas:  N„  10  cm1  min'1 

Temperature:  60°C 

Detector:  Flame  ionization 

The  retention  time  of  acetone  in  the  chromatographic  column  under  these  conditions  was  about 
3  minutes.  Although  several  peaks  in  addition  to  that  due  to  acetone  in  soil  air  were  observed  in 
the  chromatogram  at  high  sensitivity  none  of  these  peaks  interfered  with  the  acetone  determination. 

Results  and  discussion 

Adsorption  isotherms  for  acetone  adsorption  at  25°C  by  the  field-moist,  air-dry  and  oven-dry 
soil  samples  are  shown  in  Figure  1,  in  which  the  weight  of  acetone  adsorbed  per  unit  weight  of 
dry  soil  is  plotted  against  the  weight  of  acetone  per  unit  volume  in  the  vapor  phase  at  equilibrium. 
In  Figure  2,  data  are  shown  which  were  obtained  using  a  set  of  samples  composed  of  2  g  of  water 
but  no  soil.  Sectors  of  the  experimental  isotherms  obtained  for  the  four  types  of  samples  are 
drawn  in  Figure  3  for  comparison  purposes. 


g  acetone// air 
a. 


Figure  1.  Adsorption  of  acetone  by  soil  from  Laguna  Joyuda,  Puerto  Rico. 
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As  shown  in  Figure  3,  the  isotherms  obtained  for  the  adsorption  of  acetone  by  the  soil  at  field 
moisture  and  for  the  partition  of  acetone  in  water,  plotted  on  a  comparable  basis,  virtually  coincide, 
showing  that  sorption  of  acetone  by  the  soil  at  this  moisture  content  can  be  accounted  for  by  the 
partition  of  acetone  in  soil  water  rather  than  by  adsorption  on  mineral  surfaces  in  the  soil.  Re¬ 
moval  of  water  from  the  soil  by  oven-drying  enhanced  the  adsorption  of  acetone  enormously  in  com¬ 
parison  with  the  other  samples.  As  illustrated  by  the  results  obtained  with  air-dry  samples,  the 
presence  of  trace  amounts  of  water  in  the  soil  significantly  reduced  the  amount  of  acetone  adsorbed 
even  though  considerably  more  acetone  was  adsorbed  by  the  air-dry  than  by  the  field-moist  samples. 
The  effect  of  moisture  content  on  acetone  adsorption  seems  more  apparent  when  the  values  of  the 
slope  taken  from  the  linear  portion  of  the  isotherms  at  low  acetone  contents  (Fig.  1)  are  compared. 
The  numerical  values  are  about  0.11,  6.5  and  56  for  the  field-moist,  air-dry  and  oven-dry  conditions, 
respectively.  Thus  it  can  be  seen  that  oven-drving  of  the  field-moist  soil  increased  acetone 
adsorption  by  a  factor  of  more  than  500. 

These  results  are  m  general  agreement  with  those  of  Parfitt  and  Mortland1'  who  also  found  that 
acetone  is  less  strongly  adsorbed  by  clay  minerals  in  the  presence  of  adsorbed  water.  Using 
infrared  spectroscopic  techniques,  these  investigators  obtained  evidence  which  suggests  that  in 
the  absence  of  water,  acetone  is  coordinated  directly  to  the  exchangeable  cations  associated  with 
the  clay  mineral  surfaces,  but  that  in  the  presence  of  adsorbed  water,  acetone  is  weakly  bound  by 
bridging  through  water  molecules  coordinated  to  the  cations.  The  results  of  the  present  experiment 
indicate  that  at  higher  moisture  content,  the  exchangeable  cations  in  soils  are  completely  shielded 
from  acetone  by  water  molecules  and  that  the  interaction  can  be  accounted  for  by  simple  partition  of 
the  acetone  between  vapor  and  liquid. 

In  order  to  judge  the  effect  that  adsorption  of  volatile  chemicals  by  soil  might  have  on  attempts 
to  locate  a  tunnel  at  Laguna  Joyuda  by  trace  gas  detection,  consider  the  following  examples.  First, 
adsorption  by  soil  at  the  tunnel  walls  will  decrease  the  vapor  concentration  of  compound  in  the 
tunnel  air  and,  consequently,  the  amount  of  compound  which  initially  appears  in  the  atmosphere 
tlirough  ports  and  vents.  Suppose  that  acetone  vapor  in  a  tunnel  1  meter  in  diameter  is  in  equilib¬ 
rium  with  adsorbed  acetone  in  a  layer  of  soil  1  cm  thick  at  the  tunnel  wall.  It  is  assumed  that 
diffusion  of  acetone  through  air  to  the  tunnel  wall  is  fast  in  comparison  with  both  the  diffusion  of 
acetone  up  through  the  soil  profile  and  mass  movement  of  acetone  through  ports  and  vents.  Even 
though  the  total  amount  of  acetone  is  unknown,  the  distribution  of  acetone  between  the  soil  and 
air  phases  can  be  estimated.  The  result  of  this  calculation,  assuming  that  the  soil  at  the  tunnel 
wall  is  at  the  field-moisture  content,  is  that  only  12°(,  of  the  acetone  initially  present  would  remain 
in  the  tunnel  air.  At  the  air-dry  moisture  content  only  0.2^  of  the  acetone  produced  in  the  tunnel 
would  remain  in  the  air.  Thus,  the  amount  of  acetone  appearing  in  the  atmosphere  could  depend 
greatly  on  soil  parameters  such  as  moisture  content  which  influence  adsorption  of  acetone  by  soil. 

Another  consideration  is  the  effect  of  adsorption  on  the  appearance  of  volatile  chemicals  in 
the  atmosphere  by  diffusion  up  through  t he  soil  profile.  A  detailed  analysis  of  the  extent  to  which 
adsorption  will  influence  the  rate  of  diffusion  through  soil  appears  later  in  this  report.  However, 
an  estimate  of  effect  of  adsorption  on  diffusion  can  be  made  using  available  data.  An  approximate 
relationship  between  the  diffusion  coefficient  and  the  adsorption  coefficient  is  given  by  (see  eq 
7 1-75) 


0.66  <  D,, 


(1) 
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where  Dj  is  the  diffusion  coefficient  in  cm'  see"1  of  the  compound  in  soil  with  adsorption,  D„  is  the 
diffusion  coefficient  in  cm2  sec"1  of  the  compound  in  air,  <  is  the  volume  fraction  of  air  in  the  soil, 
and  K  is  the  adsorption  coefficient  expressed  as 


g  compound/cm3  soil  air 
g  compound  adsorbed/(cm3  soil  t  soil  water) 

Expressed  in  these  units,  the  numerical  values  of  the  adsorption  coefficient  of  acetone  for  the 
field-moist,  air-dry  and  oven-dry  samples  are  about  5.4  *  10"  *,  5.5  ■  10'*,  and  6  ■  10'\  respec¬ 
tively.  Since  K  is  less  than  10*2  eq  1  can  be  approximated  by 

0.66  DolzK 

n  / .  *  > 


Thus,  interaction  of  acetone  with  the  dry  soil  could  decrease  the  diffusion  coefficient  from 
that  in  air  by  a  factor  higher  than  6  >  104.  The  effect  of  air-drying  the  field-moist  soil  is  to 
decrease  the  diffusion  coefficient  of  acetone  in  soil  by  about  102.  However,  the  moisture  content 
also  influences  the  value  of  <.  Although  the  value  of  D  t  tends  to  increase  with  increasing  moisture 
content  due  to  less  adsorption,  the  rate  of  increase  in  Dy  decreases  with  moisture  content  since 
the  value  of  e  is  decreasing  as  water  enters  the  soil  pores.  In  any  event,  it  seems  clear  that 
partition  of  compounds  in  soil  water  or  adsorption  by  soil  can  significantly  reduce  the  rate  of  move¬ 
ment  upward  through  the  soil  profile. 
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Introduction 

The  diffusion  coefficient  of  a  volatile  chemical  in  soil  depends  upon  soil  water  content,  void 
porosity  and  the  degree  of  interaction  between  the  compound  and  mineral  constituents.  The  experi¬ 
ments  described  in  this  section  were  conducted  to  evaluate  the  relative  importance  of  the  above 
factors  in  determining  the  magnitude  of  the  diffusion  coefficient.  Acetone  was  selected  for  use  in 
the  study  because  it  is  a  volatile,  moderately  polar  compound  which  is  of  practical  importance  for 
personnel  detection.  In  addition,  experimental  data  are  available  on  adsorption  of  acetone  by  soil 
at  different  moisture  contents.  Soils  from  the  detection  sites  in  Puerto  Rico  were  utilized  to  provide 
background  information  applicable  for  future  work  in  the  field.  The  soils  selected  represented  a 
range  in  soil  texture  and  mineralogy.  Moist  and  water-saturated  samples  were  studied  since  these 
two  situations  represent  practical  extremes  for  the  effect  of  water  content  on  trace  gas  diffusion. 

Experimental  procedure 

A  few  of  the  pertinent  properties  of  the  six  soils  used  in  this  study  are  summarized  in  Table  I. 
The  primary  consideration  in  the  selection  of  the  soils  was  texture  which  ranged  from  sand  to  clay. 
More  detailed  descriptions  are  available  of  both  the  profile  characteristics21  and  chemical  properties22 
of  the  soils. 

Diffusion  of  acetone  through  the  soils  was  investigated  using  the  type  of  diffusion  cell  illustra¬ 
ted  in  Figure  4.  The  diffusion  system  consisted  of  a  half-cell  containing  an  acetone  source, 
separated  by  a  perforated  plate  from  another  half-cell  filled  with  soil.  Each  half-cell  was  constructed 
of  Tygon  tubing  with  an  inner  diameter  of  5  cm  and  a  height  of  10  cm.  An  air-tight  seal  was  formed 
when  the  two  half-cells  were  coupled  together. 
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Figure  4.  Dillusion-cell. 


To  determine  the  diffusion  of  acetone  vapor  in  moist 
soil,  soil  was  packed  in  a  cell  to  a  density  of  about 
1.0  g  cm'*.  The  moisture  content  of  the  soil  was  15%  of 
the  dry  weight.  The  acetone  source  consisted  of  a  solution 
of  2%  by  volume  acetone  in  water.  The  perforated  plate 
supporting  the  soil  was  separated  from  the  acetone  solution 
by  about  3  cm  so  that  the  solution  did  not  come  in  contact 
with  the  soil.  The  source  solution  was  constantly  stirred 
to  make  certain  that  there  was  no  concentration  gradient 
between  the  soil  and  the  source  solution.  Acetone  in  the 
soil  atmosphere  was  determined  by  gas  chromatographic 
analysis  (described  earlier)  of  vapor  samples  withdrawn 
from  the  cell  through  the  septums  depicted  in  Figure  4. 
Samples  were  taken  as  a  function  of  both  time  and  distance 
from  source. 

In  experiments  where  the  diffusion  of  acetone  through 
water-saturated  soil  was  investigated,  the  cells  were 
packed  with  soil  at  a  density  of  about  1.8  g  cm**.  The  soil 
in  the  cell  was  water-saturated  by  upward  percolation  of 
water.  The  water  content  of  these  samples  ranged  from  15 
to  23  g  of  water  per  100  g  of  soil  depending  upon  the  texture 
and  bulk  density.  The  same  acetone  source  as  described 
above  was  utilized.  Diffusion  was  allowed  to  proceed  into 
the  samples  for  about  24  hours.  The  soil  was  then 
microtomed  into  segments  3  to  5  mm  thick  and  stored  in  vials. 
Water  was  added  to  each  vial  to  increase  the  volume  of  the 
soil  water/acetone  solution  for  sampling  so  that  the  actual 
amount  of  acetone  in  the  soil  water  could  be  determined. 

A  24-hour  equilibration  period  was  allowed  after  addition  of 
water  prior  to  analysis  of  the  soil  solution  for  acetone  using 
gas  chromatography. 

One  additional  experiment  on  diffusion  of  acetone  in 
water-saturated  soil  was  performed  using  another  method 
which  differed  from  the  one  just  described  in  that  a  1-cm 
segment  of  soil  was  sandwiched  between  the  usual  acetone 
vapor  source  and  a  layer  of  water  in  the  diffusion  cell.  The 
water  on  top  of  the  soil  acted  as  a  trap  for  acetone  diffusing 
up  through  the  water-saturated  soil.  The  amount  of  acetone 
appearing  in  the  water  trap  was  monitored  after  a  steady-state 
condition  was  reached.  The  water  was  changed  frequently 
to  prevent  back-diffusion  of  acetone. 


Calculations 

Migration  of  volatile  chemicals  in  soil  can  be  considered  a  combined  process  of  diffusion  in  both 
the  gas  and  liquid  (water)  phases  during  which  a  relationship  exists  between  the  concentration  of  the 
chemical  in  the  gas  phase  and  liquid  phases,  and  the  amount  of  chemical  adsorbed  on  soil  surfaces. 
As  was  demonstrated  by  the  experimental  work  on  acetone  adsorption,  when  the  soil  particles  are 
covered  by  a  film  of  water  the  effect  of  adsorption  is  negligible.  This  is  the  situation  prevailing  in 
the  present  work.  We  shall  classify  our  experimental  results  into  two  categories  according  to  the 
mechanism  of  migration,  namely,  1)  moist  soil  and  2)  water-saturated  soil. 
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Dillusion  in  moist  soil.  In  this  case,  since  the  dif'fiisivity  of  acetone  in  air  is  of  a  much 
greater  magnitude  than  that  in  water,  the  rate  of  migration  is  dominated  by  diffusion  in  the  gas 
phase.  In  general,  the  relationship  between  the  partition  of  acetone  in  the  gas  and  in  the  li quid 
phases  is  in  dynamic  equilibrium.  However,  diffusion  is  thought  to  be  much  slower  Ilian  partition 
of  acetone  in  thin  films  of  water.  Therefore,  we  shall  assume  that  a  local  static  equilibrium 
condition  holds  throughout  the  moist  soil. 

The  equation  describing  diffusion  is  written  as* 


cX:;  d2c. 

*  n  * 


where 

Cj  concentration  of  compound  (g  in  the  soil  gas  phase  per  cm1  soil) 
X  distance  in  cm 

t  time  in  sec. 

The  initial  and  boundary  conditions  are: 

C,  =  0  at  t  =  0 

Cj  =  CQ  X  0  t  >  0 

C|  E  0  X  au  t  >  0  . 

The  solution  of  eq  1  with  the  conditions  eq  4-6  is  given  as 


CD 


(I) 

(5) 

(6) 


(7) 


Diffusion  in  water-saturated  soil.  The  migration  of  volatile  chemicals  is  governed  in  this  ease 
by  diffusion  in  soil  water.  Since  the  soil  pores  are  water-saturated,  the  influence  of  partition  be¬ 
tween  soil  vapor  and  water  phases  can  be  neglected.  The  diffusion  equation  is  then  written  as 


X 


d2C 
dX 2 


(8) 


where 

C  concentration  of  acetone  (g  cm*1  soil) 

D„  diffusion  coefficient  (cmJ  sec'1)  in  soil  water. 

<1 

The  solution  of  eq  8  is  the  same  as  that  of  eq  3.  Since  the  soil  contains  no  compound  at  the  be¬ 
ginning,  the  initial  and  boundary  conditions  are  written  as 


C  0 


at 


l  0 


(9) 
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A'  0 

t  0 

(10) 

c 

0 

X 

l  0. 

(11) 

Equation  10  implies  that  at  the  source  of  the  compound  the  concentration  is  held  constant.  Equation 
11  holds  if  the  experiments  are  conducted  so  that  diffusion  to  physical  boundaries  is  not  significant. 
The  solution  of  equation  S  with  the  above  initial  and  boundary  conditions  is  given  as 


where  the  error  function  is  defined  as 

x 

erf(.Y)  f  e‘u  du .  (13) 

\  n  ■' 
v  0 

For  both  the  water-saturated  soil  and  moist  soil  conditions  described  above,  the  diffusion 
coefficient  can  lie  determined  bv  comparing  experimental  results  with  those  theoretically  predicted. 
Namely,  after  obtaining  an  experimental  concentration  profile  C,.xp  (X,  t)  the  theoretical  profile 
C’Ih(.Y,  t,  D)  which  best  fits  the  experimental  profile  is  determined  by  changing  the  value  of  D. 
Although  it  is  possible  to  express  quantitatively  the  degree  of  agreement  between  two  profiles  by 
introducing  the  distance  between  two  functions  as  in  the  mathematical  analysis  of  a  functional 
space,  the  determination  of  diffusion  coefficients  in  this  work  was  made  by  visual  comparison  of 
the  experimental  and  theoretical  concentration  profiles. 

In  the  experiment  on  diffusion  of  acetone  in  water-saturated  soil  in  which  acetone  was 
monitored  in  water  above  the  soil  column  after  attaining  a  steady-state  condition,  the  diffusion 
coet ficient  from  the  experimental  data  was  calculated  using  Fick’s  law: 

'  u- 

where  the  flux  F  is  grams  of  acetone  diffusing  per  unit  of  tirea  per  second. 

Results  and  discussion 

A  comparison  of  the  experimental  plot  of  Cj  C0  versus  distance  from  the  source  after  diffusion 
of  acetone  in  moist  soil  from  Laguna  Joyuda  for  5  hours,  with  curves  obtained  utilizing  different 
values  for  D  j  in  eq  7,  is  shown  in  Figure  5a.  The  closest  agreement  was  obtained  using  a  value  of 
Dj  equal  to  3.5  •  10*\  Virtually  the  same  value  of  D,  was  also  observed  for  experimental  time 
periods  of  3.3,  1.7  and  1.1  hours.  The  experimental  and  calculated  values  for  C  j/C0  versus  distance 
for  those  results  not  shown  in  Figure  5a  are  found  in  Table  II.  These  experimental  values  for  vapoi 
diffusion  compare  very  favorably  to  t lie  one  estimated  (1.65  x  10’*  cm2  sec*1)  using  eq  1  which  indi¬ 
cates  the  validity  of  this  relationship  in  estimating  diffusion  coefficients. 

The  comparison  of  experimental  values  of  C/C0  versus  distance  (Fig.  5b)  using  water- saturated 
soil  from  Laguna  Joyuda  with  values  calculated  using  eq  12  showed  a  best  fit  using  a  value  of  Da 
of  about  1  •  10*\  Similar  values  of  D„  were  observed  for  diffusion  of  acetone  in  the  other  water- 

u 

saturated  soils  (Table  III).  This  similarity  in  values  of  the  diffusion  coefficients  obtained  using 
clav  soils  as  compared  with  those  obtained  using  soils  much  lighter  in  texture  and  lower  in  surface 
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b.  Water-saturated  soil. 


Figure  5.  Plot  ol  Cj/C0  vs  distance  from  source  for  diffusion  of  acetone  in  Laguna  Joyuda  soil. 
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Table  II.  Comparison  of  experimental  and  theoretical  values  (or 
diffusion  of  acetone  In  moist  Lacuna  Joyuda  soil. 


t  (hours) 

X  (cm) 

(C  l/(Vexp 

(C\/C0>th 

2.3 

1 

0.4 

0.62 

2.3 

2 

0.26 

0.36 

2.3 

3 

0.11 

0.14 

2.3 

4 

0 

0 

1.7 

1 

0.3 

0.36 

1.7 

2 

0.13 

0.17 

1.7 

3 

0.02 

0.02 

1.1 

1 

0.135 

0.2 

1.1 

2 

0.025 

0.03 

Table  M.  Diffusion  coefficient  of  acetone  Is  water-saturated  soils. 


Site 

Texture 

D&  cm1  sec'1 

Rio  Ouayanes 

Sandy  loam 

2  x  10‘5 

Sabana  Hoyos 

Silt  loam 

1  x  10'5 

Laguna  Joyuda 

Silt  loam 

1  x  10’5 

Lajas 

Clay 

0.9  x  10's 

Toro  Negro 

Silt  loam 

1.5  x  10*’ 

Moca 

Clay 

0.45  x  10'5 

area  supports  the  view  that  adsorption  on  mineral  surfaces  is  not  a  factor  in  influencing  diffusion 
in  moist  or  wet  soils.  Differences  in  diffusion  coefficient  are  probably  due  to  variations  in 
moisture  content  and  in  internal  geometry  of  the  samples,  the  clay  soils  having  the  lowest  diffusion 
coefficient  and  the  light-textured  soils  the  highest  value. 

Although  not  shown,  values  of  the  diffusion  coefficient  in  water-saturated  Moca  and  Lajas  soils 
calculated  using  eq  14  were  lower  than  those  reported  in  Table  III  by  only  a  factor  of  2.  The 
agreement  between  the  two  methods  lends  added  assurance  to  the  validity  of  the  values  obtained. 

The  results  of  this  study  indicate  that  vapor  diffusion  in  soils  is  not  significantly  influenced 
by  variation  in  soil  texture.  The  interaction  of  a  volatile  compound  with  moist  soil  is  primarily 
partition  in  soil  water  rather  than  adsorption  at  soil  surfaces.  The  most  important  factors  governing 
vapor  diffusion  in  soil  are  soil  water  content  and  void  porosity.  For  acetone  the  diffusion  coeffi¬ 
cient  decreased  10-  to  20-fold  upon  water-saturating  moist  soil.  This  indicates  that  appearance  of 
detectable  vapors  in  the  atmosphere  above  water-saturated  soil  will  be  considerably  slower  than  in 
the  case  of  moist  soil. 
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PREDICTION  OF  VAPOR  DIFFUSION  IN  SOIL 

Introduction 

The  experimental  work  described  in  the  previous  sections  has  been  restricted  to  laboratory 
investigation.  However,  using  these  results  it  is  now  possible  to  predict  the  influence  of  soil  on 
trace  gas  diffusion  in  the  field.  This  is  accomplished  by  first  constructing  a  model  for  calculating 
the  appearance  of  vapors  in  the  atmosphere  based  upon  one-dimensional  diffusion  from  a  subsur¬ 
face,  constant  source.  The  feasibility  of  detecting  buried  explosives  and  tunnels  is  discussed 
based  upon  this  analysis.  The  results  of  these  predictions  are  supported  by  experimental  work  on 
diffusion  of  acetone  in  the  field  at  the  Laguna  Joyuda  site  in  Puerto  Rico. 


Construction  of  model 

To  simplify  the  mathematical  analysis  for  predicting  the  influence  of  soil  on  detection  it  is 
necessary  to  introduce  assumptions  which  are  physically  acceptable.  The  process  of  primary 
concern  is  the  migration  of  volatile  chemicals  from  buried  explosives  or  a  tunnel  source  vertically 
upwards  to  the  soil/atmosphere  interface.  The  surface  of  the  source  can  be  approximated  as  a 
plane  surface  and  the  process  treated  as  one-dimensional  diffusion.  In  addition,  it  is  assumed  that 
the  soil  is  homogeneous  with  respect  to  those  properties  which  affect  the  diffusion  of  volatile 
chemicals.  The  next  concern  is  the  physical  condition  of  the  atmosphere  close  to  the  earth’s  sur¬ 
face  where  it  is  assumed  that  a  boundary  layer  exists  in  which  only  molecular  diffusion  occurs.  The 
concentration  of  compound  in  the  atmosphere  at  the  upper  surface  of  the  boundary  layer  is  taken  to 
be  negligible  relative  to  that  at  the  soil/atmosphere  interface.  A  schematic  drawing  depicting  the 
situation  described  above  is  shown  in  Figure  6.  The  ^-coordinate  is  taken  to  be  positive  perpen¬ 
dicular  to  the  earth’s  surface  in  the  direction  of  diffusion. 


ATMOSPHERE 
IMu  nq  ond  Convection) 

ATMOSPHERIC  BOUNDARY  LAYER 
I  Diffusion  Onlyl 


The  variables  not  previously  defined  which 
are  involved  in  the  analysis  are  defined  as  follows: 

s  Cj,  =  concentration  of  compound  in 

boundary  layer  air  (g  cm-1) 

S  =  concentration  of  adsorbed  compound, 
g  cm**  soil 


HOMOGENEOUS  EARTH  MEDIA 


SOURCE 

(Mme  or  Tunnel) 


L 


Ds  =  diffusion  coefficient  of  compound 
in  soil  in  absence  of  adsorption 
effect  8  (cmJ  sec*1) 

5  =  thickness  of  boundary  layer  at 
earth’s  surface  (cm) 


Figure  6.  Schematic  diagram  of  model  (or 
predicting  vapor  diffusion  from  the  surface 
of  a  buried  mine  or  a  tunnel. 


L  =  distance  between  upper  surface  of 
the  source  and  the  soil/atmosphere 
interface  (cm). 

If  local  equilibrium  can  be  assumed  to  exist  be¬ 
tween  the  gas  phase  and  immobile  (liquid  and  solid) 


phases,  the  equation  of  diffusion  of  a  volatile  chemical  through  the  earth  media  and  the  boundary 
layer  can  be  written,'  respectively, 


dCl  a2c, 

— -  =  Dj -  -  L  <  X  <  0  (15) 

dt  dX2 


14 


DETECTION  OF  BURIED  EXPLOSIVES  AND  TUNNELS 


— -  D .. 

dt  * 


0  ^  X  <  5. 


In  eq  15,  the  value  of  the  diffusion  coefficient  Dj  is  given  by 


k  *  1 


k  -  . 
C, 


The  constant  k  is  the  adsorption  coefficient  of  the  volatile  chemical  by  the  earth  material  which 
can  be  calculated  from  the  linear  portion  of  an  experimental  adsorption  isotherm  at  low  concentra¬ 
tion  levels. 

Now,  if  the  volatile  chemical  compound  under  consideration  is  not  initially  present  in  the  soil, 
then  the  initial  condition  may  be  stated  as 

C,  •  Ca  0  t  <_  0.  (19) 

If  a  volatile  compound  is  then  introduced  into  the  soil  from  the  surface  of  a  mine  or  tunnel  and  its 
concentration  at  the  source  surface  kept  constant  at  CQ  by  some  mechanism  for  t  >  0,  the  boundary 
conditions  are 


Cj  c0 


X  -  L 


C„  0 


X  =  8 


X  --  0 


X  =-  0. 


In  this  treatment,  use  of  dimensionless  variables  is  worthwhile  in  order  to  reduce  the  number 
of  variables  and,  also,  to  demonstrate  the  interdependence  of  sets  of  variables.  The  latter  point  is 
particularly  helpful  in  this  case  since  relative  values  can  be  calculated  even  though  experimental 
values  of  many  of  the  variables  are  not  yet  available. 

With  introduction  of  the  following  dimensionless  variables, 


L 


(26) 
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D* 

r  D*i 

L" 

eq  15,  16  and  19-23  become 


<)C[ 

dL’cj 

D\ 

-  1  X* 

0 

lit • 

fix'- 

<Pc!, 

— 

D* 

0  x'  ^ 

/5‘ 

tlx '  ~ 

c; 

0,  C 

'•  0 

for  l  f  0 

C1 

1 

at  xf 

-  1 

f 

0 

c f 

0 

at  x1 

t* 

:  0 

c f 

ct 

<C2 

at  x* 

0 

t1 

0 

MC 

1 

<IC A 

O'  .i 

K  — 

at  x’ 

0 

dx 

k 

ilx 1 

If. 


<2li) 


(28) 


(20) 


CUM 


CM) 


CI2) 


CUM 

CM) 


(35) 


In  eq  27  and  28,  any  arbitrary  v.alue  can  be  selected  for  D*.  a  reference  diffusion  coefficient. 
In  this  treatment,  however,  D *  cancels  in  all  final  equations  so  assignment  of  a  value  is  not  re¬ 
quired. 

Applying  the  Laplace  transform  to  eq  29-35,  one  obtains 

PC.  ^  D\  —I  (36) 

dx(i! 

(Pc' 

PC.,  D,\  --  -  (37) 

dxu> 
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at  x*  1 


<\.  0 


al  x'  .V 


C-j  «t\,  at  x* 


at  x '  0 


c  I"  c;  dt' 


U  1.2). 


The  solutions  for  the  above  system  of  equations  are  given'*  as 


eosh  </ j  x+ sinh  fl;,  <V  -  </  sinh  <?[  x4  cosh  qL,  S' 
P(cosh  sinh  <V  *  <»  sinh  ql  cosh  8‘) 

rl  sinh  t/o  (S'  -  x4) 

Pleosh  qt  sinh  q2  6'  .  o  sinh  Qj  cosh  Q2  fi4) 


1  /5S5T 


<?  i  fP/D  i  (46) 

9-j  /P/Dy  .  (47) 

Since  Cj  and  C2  are  single-valued  functions  of  the  complex  number  P  with  simple  poles,  Cj  and 
C! ,  can  be  obtained  from  eq  43  and  44  by  use  of  the  inversion  theorem.  The  result  is  that 

r<  -<S  \d  D0  -  x'  ^  ^  sin'-'  VD  fln  sin  Hn  ( 1  *x‘)exp(  /ijj/)Jt‘) 

<<VvDD0  4  1  /)„  (sin*  54  VS  4  S4  xDD^rW/y 


I'OS  (flfn  ff)  .  |  .  ,  o  <>  n  - 

- — sin  (nm it x  )  exp (  D\nim‘mnii') 


nlr  4  mo) 
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^  \1)D0(8'  -x1)  ~  sin  ?>'  yjD  (in  sin  fin  sin  XD  (/>'  x '  )/lflexp( -/ijjDjl'l 

( <V  VD  D0  <  1  Tl  *0„  (sinV  \0/i„  .  V  \V  D~' ,  1  sin"  /*„) 


cos  (nm  n)  nmrrxH 


i  n(r  ♦  mo) 


sin - t*xp  (  D'^n"  n- 1' ) 

S' 


where 


\  r 

S'  \  /  —  —  ;  r,  m  are  integers 


/fn's  are  the  positive  roots  of  the  equation 

cot  pn  i  o  cot  S'  y/UfiH  0.  (52) 

For  application  of  the  model  to  the  problem  of  detection,  it  is  only  necessary  to  use  eq  48  to 
calculate  concentration  profiles  in  the  soil,  the  concentration  in  the  air  at  the  soil/atmosphere 
interface,  and  the  flux  of  chemical  at  the  soil/atmosphere  interface.  In  order  to  calculate  the 
transient  concentration  profile  at  any  time,  the  two  infinite  series  in  eq  48  must  he  evaluated. 
Equation  48  can  be  rewritten  as 


c;  (x\  t+) 


C  J  (x\  «)  -  s  An  -  1  B„ 


where 


C\  (X  +  ,  oc)  = 


»S'  yjT)  Dq  -  X4 
<S4  y/U  Dq  i  1 


2  sin8  S4  sin /!„  (1  4  x4-)  exp  (~(i*D\P) 
(sin8  S'  \75  fin  .  8\T)  D^1  T1  sin2  (in) 


D  2  o  cos  (nm  n)  +.  ,  nt  2  2  1? 

B  -  .  -  sin  (nmrrx*)  exp(-D]nKm‘- nK  (  ). 

n(r  i  mu)  n  1 
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The  order  of  magnitude  for  An  and  Bn  has  to  be  estimated  in  order  to  decide  the  number  of 
terms  which  must  tie  computed.  From  eq  55  and  56  it  can  be  shown  that 

An  ■  2/i-‘eipl  ll;D[C)  (57) 

B  —  _  exp  (  DlirnrW).  (58) 

1  nm  n  1 

Since 

(n  -  1)  n  «  fin  riff  n  1,2...  (59) 


the  number  of  terms  to  tie  computed  for  a  given  accuracy  can  be  estimated  from 


n  1 


—  exp  |  (nJ  t)  n~  D\  t '  | 

n  l 


and 


Bn  i 

— -  -  exp  l  Ur  1  )m~  ff~Dt  tf|. 

B{  n  1 


(60) 


(61) 


In  eq  60  and  61  the  rate  of  convergence  increases  as  tf  increases.  The  roots  of  eq  52  (ftn)  can  be 
computed  as  follows  when  £*  \T)  >s  ■  mall.  Since 


cot  x 


x" 

45 


(-1)0-1  22"  B 


(2n) 


x2n-l 


(62) 


where  Hn  is  the  Bernoulli  number,  it  follows  that 

cot  y  VD  Pn  («*  »f  ^  V^/8n  <•  1  (63) 


and  that 


cot  fln 


n  rr 


Thus,  from  eq  52 


n 


n  rr 

D  t-  1 


(64) 


(65) 


If  8*  \Z5  is  not  small  enough,  the  Newton-Raphson’s  method  can  be  used  to  obtain  the  roots. 

From  eq  48  it  is  seen  that  the  steady-state  concentration  profile  at  infinite  time,  Cj  (x*\  *>), 
is  given  by  eq  54.  Thus,  the  maximum  obtainable  concentration  at  the  soil/atmosphere  interface 
at  x*  0  can  be  calculated  from 
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(8*  yjD  D0 

q(0,  ~)  =  - .  (66) 

(8^^/DDq  +  1 

The  relative  concentration  Cf  indicating  the  degree  of  approach  to  the  steady-state  condition  at 
t  =  »  is  given  by  the  ratio  of  eq  48  to  eq  54. 

C t  (x+,  t^)  Cj(*+,  t+) 

Cf(x\  t4)  - -  -  .  (67) 

C;(x+.  «)  Cjfx4.-) 

The  flux  of  chemical  at  the  soil/atmosphere  interface  is  given  bv 


where  F  is  the  flux  in  g  cm'1  sec*'.  Using  the  dimensionless  parameters  previously  defined,  one 
obtains 


The  concentration  gradient  (dC \/dx+)  can  be  evaluated  using  eq  48. 

Computations 

Using  the  model  developed  in  the  previous  section,  it  is  now  possible  to  predict  the  rate  of 
movement  of  trace  gases  originating  from  mines  and  tunnels  and  moving  up  through  soils.  As  an 
example,  the  movement  of  acetone  through  soil  will  be  examined  by  incorporating  into  the  model 
results  of  investigations  of  the  diffusion  of  acetone  through  soil  from  the  detection  site  at  Laguna 
Joyuda,  Puerto  Rico.  Inasmuch  as  acetone  is  given  off  by  people,  the  results  obtained  should  be 
directly  applicable  to  detection  of  occupied  tunnels.  Although  acetone  is  probably  not  associated 
with  mines,  TNT  vapors  are  expected  to  interact  less  than  acetone  with  moist  soils,  as  will  be 
explained  later.  Thus,  acetone  diffusion  from  a  source  located  near  the  soil  surface  may  be  slower 
than  the  diffusion  of  TNT  vapors.  In  any  event,  these  calculations  should  provide  a  reasonable 
indication  of  whether  mines  and  tunnels  can  be  located  by  detection  of  trace  gases  which  have 
diffused  through  soil. 

First,  the  range  of  variables  appearing  in  eq  48  should  be  established.  The  value  of  <  in  a  moist 
soil  is  likely  to  be  in  the  range 

0.1  <  '  <0.5  (70) 

with  the  exact  value  depending  on  soil  texture  and  water  content.  The  relationship  between  the 
diffusion  coefficient  Dj  and  the  adsorption  coefficient  k  is  given  in  eq  17.  The  relationship 
between  the  value  of  k  and  that  of  K,  the  adsorption  coefficient  previously  discussed  (eq  1),  is 


k 


(1  -  0 

(K 


(71) 


r 
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Now,  le'  us  assume  that  the  minimum  value  of  K  expected  is  that  which  corresponds  to  maximum 
adsorption  of  acetone  by  air-dry  Laguna  Joyuda  soil.  If  the  maximum  value  for  K  indicates  no 
adsorption,  then 


5.5  -  10~5  •_  K  (72) 

The  value  of  09  for  acetone  in  eq  17  can  be  estimated  using  the  relationship1 

Df.  0.66  <D2  (73) 

where  as  a  general  rule14 

D2  '  0.1  cm1  sec*1.  (74) 

Combining  eq  17,  71  and  73, 

0.66 

D.  -  - — — - -  .  (75) 


Equation  75  is  the  same  as  eq  1  which  was  shown  through  the  experimental  work  described  earlier 
to  be  valid  for  calculating  the  diffusion  coefficient  of  vapor  in  soil.  Thus,  considering  conditions 
70  and  72: 

5.1  x  10-8  <  Dt  1  0.033.  (76) 

The  range  in  S+  is  difficult  to  estimate  since  the  thickness  of  the  boundary  layer  8  is  unknown, 
but  if  the  depth  of  a  source  ranges  from  10  to  100  cm  and  the  thickness  of  the  boundary  layer  is 
1  cm,  then 

10~2  <  <0.1.  (77) 

Examination  of  eq  48  shows  that  C  J  is  a  function  of  xf  and  Djl+  at  constant  values  of  t,  k 
or  K,  and  Values  of  Cj  (x+,  t+)  have  been  calculated  for  the  extremes  in  t,  K  and  8+  indicated 
above  in  conditions  70,  72  and  77  and  are  tabulated  in  Tables  IV  and  V.  In  addition,  the  correspond¬ 
ing  values  of  Cf(x\  t+)  as  defined  in  eq  67  are  shown.  The  concentration  gradient  in  eq  69  for  the 
transient  condition  at  the  soil/atmosphere  interface,  G(0,  t+),  and  the  relative  concentration  gradient 
indicating  the  degree  of  approach  to  the  steady-state  condition,  G v  are  also  given  in  Tables  IV  and 
V  as,  respectively. 


G(0,  t+) 


G(0,  t+) 


(78) 


G(0,  *>) 


(79) 


DETECTION  OF  BURIED  EXPLOSIVES  AND  TUNNELS 


21 


Table  IV.  Values  of  C|(x+,  t+)  and  ',,(x+,  t  *")  predicted  for  extreme 
(low)  values  of  r ,  S+,  sad  K. 

€  -=  0.1;  S4  -  10TZ ;  K  5.5  x  kT5 


0{t4  =  0.04  Djt+  =  0.08  D\t'  0.18 


x+ 

Cj(x+,  t+) 

Cjfx4’,  t+) 

Cj(x+,  t+) 

Cjlx' ,  t4) 

C\(x\  l4i 

C/x\  I4) 

-0.9 

0.724 

0.804 

0.803 

0.892 

0.838 

0.931 

-0.8 

0.480 

0.599 

0.617 

0.771 

0.683 

0.854 

-0.7 

0.289 

0.413 

0.453 

0.648 

0.540 

0.771 

-0.6 

0.157 

0.262 

0,317 

0.529 

0.414 

0.689 

-0.5 

0.077 

0.154 

0.211 

0.422 

0.305 

0.610 

-0.4 

0.034 

0.085 

0.133 

0.333 

0.216 

0.541 

-0.3 

0.013 

0.044 

0.079 

0.263 

0.145 

0.484 

-0.2 

0.005 

0.023 

0.043 

0.214 

0.088 

0.441 

-0.1 

0.001 

0.014 

0.018 

0.185 

0.042 

0.414 

0 

1.74  x  Kf* 

0.011 

2.74  x  10-8 

0.178 

6.27  x  10-8 

0.407 

0(0,  t) 

0.011 

0.175 

0.406 

°r 

0.011 

0.175 

0.406 

£>{t+  =  0.16 

Djt+  =  0.20 

Ojl4  - 

-0.9 

0.869 

0.955 

0.873 

0.970 

0.900 

1.000 

-0.8 

0.722 

0.903 

0.748 

0.935 

0.800 

1.000 

-0.7 

0.593 

0.948 

0.628 

0.900 

0.700 

1.000 

-0.8 

0.475 

0.791 

0.676 

0.860 

0.600 

1.000 

-0.5 

0.369 

0.738 

0.412 

0.823 

0.500 

1.000 

-0.4 

0.276 

0.689 

0.316 

0.790 

0.400 

1.000 

-0.3 

0.194 

0.648 

0.229 

0.762 

0.300 

1.000 

-0.2 

0.123 

0.617 

0.148 

0.741 

0.200 

1.000 

-0.1 

0.060 

0.600 

0.073 

0.727 

0.100 

1.000 

0 

9.11  x  10-* 

0.593 

1.11  x  10"7 

0.724 

1.54  x  10-7 

1.000 

0(0.  t) 

0.591 

0.723 

1.000 

°r 

0.691 

0.723 

1.000 

Upon  comparing  the  calculated  values  in  Tables  IV  and  V,  it  is  clear  that  r ,  K  and  <5  4  in  the 
range  of  interest  have  little  effect  on  the  values  of  C \  and  C  for  x+  <  -0.1.  This  is  because  for  a 
given  value  of  Djt+,  x+  becomes  the  predominant  variable.  For  example,  for  (<,  K,  8*)  <  (0.5,  «, 
0.1),  in  eq  54  the  term  {8+  y/D  DQ  <  1.65  x  10-2.  Therefore,  eq  53  can  be  approximated  by 

C1  =  *i  ~  %  An  ~  *  Bd  .  (80) 

n=l  n=l 

As  the  value  of  x+  approaches  zero,  the  effect  of  the  other  variables  on  C  J  begins  to  predominate. 
However,  the  values  of  Cf(x+,  t+)  and  G(0,  t+)  are  insensitive  to  variation  in  <,  K,  and  S4  even  at 
x+  =  0.  Inasmuch  as  the  values  in  Tables  IV  and  V  are  applicable  over  the  entire  range  of  the 
variables  of  interest,  they  have  been  summarized  in  Figures  7-9  in  which  Cj(x\  f4)  and  Cr(x+,  t+) 
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Table  V.  Values  of  C{(x4,  r4)  and  Cr(x4,  t4)  predicted  for  extreme 
(high)  values  of « ,  <  ,  and  K . 


X* 

i  =  0.5;  0.1 

,  K  =  * 

DjV 

0.04 

fl[t4  =  0.08 

1 

•-  + 

-t- 

II  1 

0.12 

Cj(x\  t') 

c/x\  C) 

C4(x4,  t4) 

Cjfx4,  tfl 

C{fx\  t4! 

C/x4.  t+) 

0.9 

0,724 

0.800 

0.803 

0.890 

0.838 

0.930 

0.8 

0.480 

0.597 

0.617 

0.768 

0.683 

0.850 

-0.7 

0.289 

0.410 

0.453 

0.043 

0.540 

0.700 

-0.0 

0.157 

0.259 

0.317 

0.523 

0.413 

0.082 

-0.5 

0.077 

0. 152 

0.211 

0.410 

0.308 

0.602 

-0.4 

0.034 

0.083 

0.133 

0.325 

0.217 

0.530 

-0.3 

0.013 

0.043 

0.079 

0.254 

0.146 

0.470 

-0.3 

0.005 

0.022 

0.043 

0.204 

0.090 

0.425 

-0.1 

0.001 

0.012 

0.020 

0.173 

0.044 

0.396 

0 

1.30  x  10~4 

8.00  *  10~3 

2.53  x  10-3 

0.156 

6. 16  x  10~3 

0.380 

G0 

9.21  x 

10  3 

0.159 

0.379 

Gr 

9.30  x 

10'3 

0.102 

0.386 

D^tf 

0.16 

D|  t4  =  0.20 

D\S  = 

-0.9 

0.859 

0.953 

0.873 

0.968 

0.902 

1.000 

-0.8 

0.723 

0.900 

0.748 

0.932 

0.803 

1.000 

-0.7 

0.594 

0.842 

0.629 

0.893 

0.705 

1.000 

-0.0 

0.470 

0.784 

0.517 

0.853 

0.006 

1.000 

-0.5 

0.370 

0.728 

0.414 

0.814 

0.508 

1.000 

-0.4 

0.278 

0.077 

0.319 

0.779 

0.410 

1.000 

-0.3 

0.197 

0.634 

0.233 

0.749 

0.311 

1.000 

-0.2 

0.128 

0.601 

0.155 

0.726 

0.213 

1.000 

-0.1 

0.066 

0.579 

0.082 

0.711 

0.115 

1.000 

0 

9.20  x  10-3 

0.567 

11.4  x  ur4 

0.702 

0.102 

1.000 

G0 

0.562 

0.694 

0.984 

Gr 

0.571 

0.705 

1.000 

are  plotted  as  a  function  of  x4  for  various  values  of  Djt4,  and  G(0,  t4)  is  plotted  as  a  function  of 
(Djt4)-1,  respectively.  As  will  be  illustrated  in  the  examples  which  follow,  use  of  these  graphs 
allows  rapid  estimation  of  Ct(x4,  t4)  and  F(0,  l4)  when  t,  8t,  L,  K  and  CQ  are  known  parameters. 

Based  on  the  results  obtained  for  acetone  the  possibility  of  detecting  mine  and  tunnel  vapors 
diffusing  up  through  the  soil  at  the  proposed  mine/tunnel  detection  research  site  at  Laguna  Joyuda 
will  now  be  examined.  For  the  case  of  a  mine,  suppose  that  t  =  0.5,  L  =  10  cm,  8  =  1cm, 

C0  =  3  x  10-11  g  cm*’,  and  K  =  5.5  x  10-5,  5.4  x  10~3andoc.  The  value  of  <  was  taken  to 
be  0.5  since  the  soil  is  likely  to  have  a  relatively  high  porosity  due  to  disturbance  upon  emplace¬ 
ment  of  the  mine.  The  values  of  the  adsorption  coefficient  correspond  to  those  determined  for 
adsorption  of  acetone  by  dry  and  moist  Laguna  Joyuda  soil,  and  no  adsorption,  respectively.  A 
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Figure  9,  Concentration  gradient  at  soil/ atmosphere  interface  as  a 

function  of  time. 


depth  of  10  cm  was  selected  but  the  distance  between  the  soil  and  mine  surfaces  is  likely  to  be  less. 
The  value  of  C0  is  that  obtained  for  the  amount  of  TNT  in  vapor  in  equilibrium  with  pure  TNT  as¬ 
suming  a  partial  pressure  of  10'*  mm.1  17 

For  the  case  of  a  tunnel,  suppose  that  f  =  0.3,  L  =  100  cm  and  C0  =  1.7  x  10*7  g  cm*’,  the 
remaining  parameters  having  the  same  values  as  given  for  the  case  of  a  mine.  The  value  of  <  was 
selected  because  the  undisturbed  soil  between  the  tunnel  and  soil  surface  should  be  of  moderate 
porosity.  A  depth  of  100  cm  is  conservative  but,  as  will  be  shown,  diffusion  to  the  surface  from 
even  this  depth  is  probably  too  slow  for  detection  of  these  vapors  to  be  significant  in  tunnel  detection. 
Inasmuch  as  there  was  no  way  to  obtain  a  more  reliable  value  of  C0,  it  was  estimated  by  assuming  a 
steady  concentration  equal  to  the  amount  of  acetone  produced  by  one  person  (0.1  mg  day'*) 1  in  a  cube 
200  cm  on  edge  in  one  day. 

The  real  time  corresponding  to  the  various  values  of  D\t*  in  Figures  7-9  is  calculated  by 
combining  eq  17  ,  27,  28  and  73  to  obtain 


t 


(1  -  «) 
eK 


1  L 


0.66  « D, 


--  D{t+ 


(81) 


It  should  be  noted  that  for  most  practical  cases  where  K  <0.1, 
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t  ~  i1 - 0  f-2  D\t'  (82) 

0.66  <2D2K 

but  that  it  K  ±  10. 

I  2 

t  ~  - - -  .  (83) 

0.66  ( D 2 

The  value  of  Cj(x+,  t+)  for  any  depth  greater  than  x^  <  -0.1  can  be  obtained  directly  using 
Figure  7  since 


Cj(x+,  t+)  =  C0Cj(x+,  t+).  (84) 

The  value  of  Cf(x+,  t+)  for  x+  <  0  can  be  taken  directly  from  Figure  8.  The  following  procedure 
can  be  used  to  calculate  the  concentration  of  compound  at  the  soil/atmosphere  interface,  C  j(0,  t). 
By  combining  eq  26  and  67  one  obtains 

Cj(0,  t+)  =  CQCj  (0,  Cr(0,  t+).  (85) 

Cf(0,  t+)  is  shown  in  Figure  8  and  Cj  (0,  «)  can  be  calculated  using 

c+  (0>  =  (86) 

1  Z  +  1 


where 


Z 


0.66  SV 


i  +  (-Lj_l)|  /j 

tK 


(87) 


Equation  86  is  obtained  by  combining  eq  17  ,  50,  66,  71  and  73.  For  most  practical  cases  where 
K  <  0.1, 


C\  (0,  ~)  = 


0.66  A  ft‘» 

(1  -  o'/j 


If  K  >  10,  then 


Cj  (0,  oc)  =  0.66  5+  <2. 


The  flux  of  compound  at  the  soil/atmosphere  interface,  F(0,  f f),  is  given  by 


F(0,  tf) 


0.66  (  D 2  Cq 


G  (0,  tf) 


(88) 


(89) 


L 


(90) 
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which  is  obtained  from  eq  69  and  78.  The  value  of  the  concentration  gradient,  G(0,  tf),  is  taken 
directly  from  Figure  9. 

Results  and  discussion 

The  calculated  values  of  C\  (0,  tf)  and  Cf  (0,  t+)  for  the  soil/atmosphere  interface  are  given 
in  Table  VI.  For  the  conditions  corresponding  to  the  case  of  a  mine  10  ctn  deep,  the  steady-state 
value  at  infinite  time  of  C  j  (0,  tf)  ranges  from  about  10'4  to  10"1  for  dry  soil  and  nonadsorbing  soil, 
respectively.  Thus,  even  if  the  chemical  compound  is  not  adsorbed  by  the  soil,  the  maximum  con¬ 
centration  expected  at  the  soil/atmosphere  interface  is  only  about  1%  of  that  at  the  mine  surface. 
For  moist  soil,  the  most  realistic  situation,  the  maximum  attainable  concentration  is  about  0.01% 
of  that  at  the  mine  surface.  Judging  from  the  values  of  Cf  (0,  t+),  it  would  require  77,  0.78,  and 
U.004  days  to  attain  a  concentration  equivalent  to  about  40%  of  the  steady-state  concentration  for 
dry,  moist  and  ideal  soil  with  no  adsorption,  respectively.  Clearly,  adsorption  of  trace  gases  by 
soil  is  an  important  consideration  in  trace  gas  detection  in  that  adsorption  markedly  decreases  both 
the  maximum  concentr  ation  possible  at  the  soil/ atmosphere  interface  and  the  rate  of  attainment  of 
the  steady-state  concentration. 

The  primary  difference  in  the  assumed  conditions  between  the  cases  of  a  mine  and  a  tunnel 
is  the  depth  L  below  the  soil  surface.  As  with  the  effect  of  adsorption,  increasing  the  depth  of  the 
source  from  10  to  100  cm  causes  both  a  decrease  in  the  value  of  Cj  (0,  *>)  and  an  increase  in  the 
time  required  to  attain  a  given  value  of  Cr(0,  l).  For  example,  considering  moist  soil,  the  ratio  of 
the  soil/atmosphere  interface  concentration  to  that  in  the  tunnel  atmosphere  at  the  steady  state  is 
predicted  to  be  about  10* 5 ,  a  factor  of  100  lower  than  that  predicted  for  the  case  of  a  mine.  Also, 
it  would  require  about  300  days  to  reach  40%  of  the  steady-state  concentration  as  compared  to  a 
time  of  less  than  one  day. 

Table  VI.  Relative  effluent  concentration  at  the  soil/atoosphere  interface  as 
predicted  for  mine  and  tunnel  installations. 
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It  seems  worthwhile  emphasizing  that  in  comparing  the  cases  of  a  mine  and  a  tunnel,  it  has 
been  assumed  that  the  thickness  of  the  boundary  layer  at  the  soil/atmosphere  interface,  8,  is 
constant.  Thus,  in  changing  the  depth  from  10  to  100  cm,  the  effect  is  to  decrease  the  value  of 
from  0.1  to  10'J.  The  importance  of  the  atmospheric  boundary  layer  thickness  becomes  apparent 
when  it  is  recognized  that  the  same  results  would  be  obtained  for  a  mine  10  cm  deep  and  a  tunnel 
100  cm  deep  if  the  corresponding  boundary  layer  thicknesses  were  1  cm  and  10  cm,  respectively. 
Similarly,  the  results  given  in  Table  VI  for  the  case  of  a  tunnel  would  be  valid  for  the  case  of  a 
mine  if  the  mine  were  10  cm  deep  and  the  boundary  layer  thickness  0.1  cm. 

The  dependence  of  C|  (0,  «.)  on  the  parameters  t,  S'  and  K  is  given  for  practical  cases  by 
eq  87.  Although  the  value  of  C|  (0,  <*)  is  most  sensitive  to  changes  in  <  followed  by  8*  and  K, 
the  range  in  these  parameters  given  by  eq  70,  72  and  77  indicates  that  in  practice,  K  and  5 4  are 
the  most  important  parameters  in  determining  the  steady-state  concentration  at  the  soil/atmosphere 
interface. 

In  Table  IV,  the  predicted  concentration  of  volatile  chemical  in  the  gas  at  the  soil/atmosphere 
interface  is  given.  The  choice  of  concentration  at  the  source,  CQ,  was  discussed  earlier.  The  flux 
of  compound  at  the  soil  surface  has  also  been  computed.  Finally,  the  time  required  to  collect  a 
minimum  detectable  amount  of  the  compound  under  ideal  conditions  has  been  estimated  in  the 
following  way.  Let 


F (0,  t4)  —  (91) 

At 

where  u>  is  the  minimum  detectable  weight  of  compound  in  the  sample  in  g,  A  is  the  soil  area 
directly  above  the  source  from  which  the  collection  is  made  in  cmJ,  and  f  is  the  sample  collection 
time  in  seconds.  Thus,  it  is  assumed  that  the  sample  is  collected  directly  above  the  source  for  a 
time  T  with  100%  efficiency.  The  value  of  w,  10'10  g,  was  selected  because  in  our  laboratory  it  is 
possible  to  detect  this  amount  of  compound  in  a  total  collection  volume  of  3  liters  of  gas  which 
corresponds  to  about  10*“  to  10’ 14  g  cm'5.  A  collection  area  of  100  and  1000  cmJ  was  selected  for 
the  cases  of  a  mine  and  a  tunnel,  respectively.  The  value  of  F( 0,  t)  is  that  calculated  using  eq  89. 

If  the  adsorption  behavior  of  TNT  were  similar  to  that  of  acetone  the  concentration  of  TNT 
in  the  vapor  near  the  soil  surface,  as  shown  in  Table  VII,  would  be  about  10*“,  10"14  and  10'“  g  cm‘s 
for  times  of  77,  0.78  and  0.004  days  after  emplacement  of  a  TNT  mine  in  dry,  moist  and  non¬ 
adsorbing  soil,  respectively.  For  all  three  conditions,  the  flux  of  TNT  vapor  from  the  soil  above 
the  mine  into  the  atmosphere  would  be  about  10'14  g  cm'J  sec'1.  The  estimated  time  necessary  to 
collect  a  minimum  detectable  sample  should  be  less  than  one  minute.  From  the  work  on  acetone 
adsorption  it  is  known  that  the  interaction  of  acetone  with  moist  soil  is  due  primarily  to  the  partition 
of  acetone  in  films  of  water  rather  than  to  adsorption  on  mineral  surfaces.  Since  TNT  is  less  soluble 
in  water  than  acetone,  one  would  expect  that  it  would  interact  less  than  acetone  with  soil  under  the 
same  moisture  conditions.  This  means  that  the  possibility  of  detecting  mines  in  the  moist  soil  by 
trace  gas  analysis  may  be  more  favorable  than  is  indicated  in  Table  VII.  These  results  support  the 
view  that  locating  mines  by  detecting  explosive  vapors  which  diffuse  through  soil  is  possible. 

For  the  case  of  a  tunnel,  the  concentration  of  compound  at  the  soil/atmosphere  interface  of  dry, 
moist  and  non-adsorbing  soils  should  be  about  10* ",  10'“,  and  10'15  g  cm'5  of  vapor  about  2  x  10\ 
200  and  0.5  days,  respectively,  after  initial  release  of  a  compound  within  the  tunnel.  The 
corresponding  flux  of  compound  into  the  atmosphere  would  be  about  10* 14  g  cm*J  sec'1.  A  minimum 
detectable  sample  might  be  collected  in  about  20  seconds.  These  results  show  that  detection  of 
trace  gases  which  have  diffused  up  through  soil  is  not  likely  to  be  of  importance  in  locating  tunnels 
by  trace  gas  analysis.  This  is  because  of  prohibitive  time  which  must  elapse  after  release  of  a 
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Table  VD.  Effluent  concentration  and  flux  at  the  toil/ atmosphere  interface,  and  minimum 
sampling  time  as  predicted  for  mine  and  tunnel  Installations. 
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compound  within  a  tunnel  before  a  detectable  atmospheric  sample  can  be  collected.  Thus,  the 
movement  of  trace  gas  through  entrances  and  ducts  should  be  of  primary  concern  in  tunnel  detection. 
However,  consideration  should  be  given  to  the  effect  of  the  interaction  of  trace  gases  with  soil  at 
the  tunnel  walls  on  the  concentration  of  effluent  which  appears  in  the  atmosphere. 

Although  no  comprehensive  experimental  investigation  of  diffusion  under  field  conditions  has 
been  undertaken,  one  preliminary  experiment  has  been  completed  which  supports  the  results  of  the 
laboratory  and  theoretical  work.  In  performing  this  experiment  10-ml  quantities  of  liquid  acetone 
were  deposited  at  depths  of  1,  2,  4  and  8  feet  in  the  soil  near  the  Laguna  Joyuda  site.  The  acetone 
was  then  allowed  to  diffuse  for  a  period  of  3  weeks  at  which  time  the  atmosphere  above  the  soil  at 
each  location  was  sampled,  and  the  acetone  concentrations  determined  by  gas  chromatography.  The 
results  (Table  VIII)  indicate  that  the  acetone  concentration  found  at  the  2-  and  4-ft  locations  was 
significantly  greater  than  that  normally  present  in  the  ambient  atmosphere.  The  data  suggest  that 
the  acetone  in  the  soil  at  the  1-ft  depth  had  already  dissipated  whereas  that  at  the  8-ft  depth  had 
not  yet  appeared  in  the  atmosphere.  These  results  support  the  view  that  trace  gas  diffusion  through 
soil  is  rapid  enough  to  be  a  consideration  in  mine  detection  but  too  slow  to  be  of  importance  in 
tunnel  detection  by  trace  gas  analysis. 

In  addition  to  acetone,  the  gas  chromatographic  analysis  indicated  the  presence  of  several 
other  compounds  present  in  concentrations  greater  than  normally  expected  in  the  absence  of  implaced 
acetone.  It  is  possible  that  these  compounds  resulted  from  biodegradation  of  the  acetone  deposited 
in  the  soil.  One  of  the  compounds  has  been  tentatively  identified  by  its  gas  chromatographic  re¬ 
tention  time  as  propane  which  could  be  formed  by  the  reduction  of  the  ketone  group  of  acetone.  The 
other  two  compounds  have  not  been  identified  but  have  retention  times  close  to  those  of  acetaldehyde 
and  2-propanol,  both  of  which  are  possible  reduction  products  of  acetone.  Although  no  quantitative 
results  are  available,  the  relative  concentrations  determined  using  chromatographic  peak  areas  are 
given  in  Table  IX. 
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Table  VIII.  Concentration  of  acetone  in  atmosphere  above  soil  containing 
acetone  emplaced  at  various  depths  at  Laguna  Joyuda. 


Depth 

Acetone  cone 

(It) 

(pg/ir 

1 

1.5 

8 

21 

4 

7.2 

8 

0.8 

No  acetone  emplaced 

0.8 

Hanover,  N.H.  air 

0.7 

*  Assumes  100%  recovery.  Checks  of  the  re¬ 
covery  efficiency  of  acetone  injected  into 
evacuated  sample  cylinders  indicate  that  about 
70%  is  recovered  by  this  method. 


Table  IX.  Relative  amounts  of  possible  decomposition  products  of  acetone 
emplaced  at  various  depths  in  soil  at  Laguna  Joyuda. 
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L0 
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4.1 

1.0 

0.0 

*Not  measurable  in  absence  of  acetone  emplacement;  Hanover  air  used  as 
reference  sample. 


A  STATISTICAL  METHOD  FOR  ANALYSIS  OF  DIFFUSION  THROUGH  SOIL 
Introduction 

When  the  mechanism  of  transport  of  volatile  chemicals  in  soils  is  governed  by  Ficlt’s  law  of 
diffusion,  theoretical  analysis  involves  solving  a  certain  differential  equation  of  a  parabolic  type 
under  properly  imposed  initial  and  boundary  conditions.  Unfortunately,  an  exact  mathematical 
solution  is  not  always  possible  for  many  practical  situations.  This  is  because  the  diffusion 
coefficient  is  seldom  constant  but  varies  with  location,  concentration  and  time  due  to  the  structural 
inhomogeneity  of  the  soil  and  to  variable  physical  conditions  such  as  temperature  and  moisture 
content  within  the  soil.  For  such  a  problem  of  variable  coefficient,  the  finite  difference  method' 
has  been  considered  one  of  the  most  effective  methods  of  analysis.  However,  there  is  another  quite 
different  approach,  the  so-called  method  of  statistical  trials  or  Monte  Carlo  method.  Using  the 
Monte  Carlo  method,  many  different  types  of  problems  of  computational  mathematics  can  be  solved  by 
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defining  each  particular  problem  in  terms  of  a  random  process.20  Boundary  value  and  initial  value 
problems  for  linear  differential  equations  are  among  the  most  interesting  fields  of  application  of 
the  Monte  Carlo  method.  The  general  mathematical  scheme  of  the  Monte  Carlo  method  for  solution 
of  differential  equations  involves  construction  of  random  walk  processes  with  finite  sets  of  states. 
For  cases  involving  diffusion  in  which  the  underlying  phenomenon  is  a  random  process,  problems 
can  be  solved  by  direct  construction  of  the  corresponding  random  process  rather  than  by  solution  of 
a  differential  equation  as  is  normally  done  by  the  Monte  Carlo  method. 

The  connection  between  the  linear  parabolic  type  of  partial  differential  equation  and  random 
processes,  so-called  stochastic  diffusion  processes,  has  long  been  known."  The  first  attempt  to 
apply  this  relation  to  the  solution  of  the  equations  was  made  by  Sragovich.21  In  the  present  work, 
a  method  based  on  stochastic  diffusion  processes  for  analysis  of  diffusion  problems  is  developed. 
The  method  is  then  demonstrated  for  acetone  diffusion  using  the  data  obtained  on  the  soil  at  the 
Laguna  Joyuda  detection  site. 

Theory 

Consider  the  following  initial-value  problem  in  three-dimensional  space: 


d C(x,  £) 

7n 


i  3 

M=1 


<)xi  dXj 


[bit(x,  t)  C  (x,  t)l  -  [a^x,  t)C  (x,  01 


(92) 


j=l 


with  the  initial  condition 
C(x,  0)  =  C0(x) 


(93) 


where  x  denotes  ajx>int  (xJf  x2,  x3).  The  problem  can  be  described20  22  in  terms  of  a  corresponding 
diffusion  process  £t(£lt,  £2t,  f3t),  with  a  diffusion  matrix  B  =  ||5^||  and  a  transition  vector  a 
(at,  a2,  a3).  The  probabilistic  meaning  of  the  solution  C(x,  t)  of  eq  92  is  the  unconditional  prob¬ 
ability  density  of  the  diffusion  process  at  time  t.  Therefore,  the  problem  reduces  to  one  of  finding 
the  probability  density  of  £v 

The  diffusion  process  can  be  constructed  from  a  Brownian  process  or  a  so-called  Wiener 
process,  qt(»?u,  q2t,  q3t),  whose  mathematical  definition  is  a  normal,  real  process  with  an  independ¬ 
ent  increment,  in  which 


E(qt  -  >7S)  =  0  •  (94) 

E(jt  -  qs)2  =  It  -  si  (95) 

where  E  denotes  a  mathematical  expectation.  Specifically,  the  relation  between  and  r/t  is 
described  in  the  following  manner  by  Ito:22  25 


(u  =  (i 


i0 


/  aj«.  fs)ds 
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1 

/=! 


r a (*• 


(96) 


where  rju,  q2t  and  q3t  are  independent  in  the  set  of  Wiener  processes,  the  function  r^(t,  £)  forms  a 
matrix  R  which  is  connected  with  B  by  the  relation 


R  •  R  B 


(97) 
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or,  in  more  detail, 

3  _ 

b^a,  a  1  rie(t,  a  rfe  ( t ,  a  1  1  i,  I  1  3  (98) 

e=l 


The  system  of  eq  96  has  a  solution  which  is  both  continuous  with  a  probability  of  one  and  unique 
under  the  condition  that  the  functions  and  satisfy  a  Lipschitz  condition.  Although  the 
Lipschitz  condition  requires  more  than  continuity,  it  does  not  impose  a  serious  restriction  on 
practical  applications.  Under  this  condition  eq  92  and  96  are  merely  two  different  mathematical 
descriptions  of  the  same  problem.  Therefore,  we  seek  the  solution  of  eq  96  instead  of  eq  92.  The 
solution  of  eq  96  may  be  constructed  by  means  of  a  direct  finite  difference  approximation  with 
respect  to  time.  In  order  to  perform  this  procedure,  we  dissect  the  interval  (0,  T)  by  the  points 

0  =  t0  <  tj  <  t2  <  ...  <  tn_t  <  t„  =  T  (99) 

into  parts  of  length 

At,  (100) 

First,  we  determine  an  initial  value,  £0  (£10,  £20,  £30),  for  the  process  from  the  given  initial 
condition,  eq  93.  The  initial  curve  CQ(x)  is  interpreted  as  the  density  distribution  of  £0.  Thus, 
the  initial  value  of  is  determined  by  the  method  of  statistical  trial  according  to  the  density 
distribution.  The  diffusion  process  £(  is  computed  by  means  of  the  random  variables 

A? j  =  1,  2,  3;  m  =  0,  1 . n  -  1  (101) 

which  are  independent  and  normal,  with  parameter  (0,  v^Atm)  from  the  following  formula: 

3 

+  ,  =  V  +  >Atm  +  2  ri.(tm,  )  A T).(m) 

m  +  l  m  m  */  ™  I 

m  =  0.  1 . n  -  1.  (102) 

With  sufficiently  fine  dissection  of  the  interval  (0,  T),  the  sequence  fta(0  <  a  <  n)  approximates 
the  continuous  process  within  any  specified  accuracy. 

There  are  several  methods  of  generating  the  random  variable  A  17  on  a  digital  computer.”  In  the 
present  work,  random  numbers  were  generated  by  a  multiplicative  congruential  procedure  in  which  the 
following  recurrence  relation  was  used: 

Rn+l  =  aRn(modm).  (103) 

From  the  many  possible  alternatives  of  the  proper  choice  of  parameters  a  and  m,  a  suggestion  by 
Moshman”  was  adapted  in  which 

a  =  2(p/S)  +  1  +  3  (104) 

m  =  yP 


(105) 
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for  a  machine  with  radix  y  having  a  word  length  of  p  digits.  The  symbol  X  denotes  the  largest 
integer  less  than  X.  The  random  real  numbers  Sa  uniformly  distributed  between  0  and  1.0  can  be 
generated  from  Rn  as 


Rn/R 
n  max 


(106) 


where  RmSLX  denotes  the  maximum  value  of  Rn.  Finally,  the  random  variable  A  q(m)  is  constructed* 
from  Sn  by  the  use  of  the  relation 


A  rj(m) 
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v 


i  1 


K 

2 


\  !(/ 12 


(107) 


in  which  A  q(m)  approaches  a  true  normal  distribution  asymptotically  as  K  increases. 


Application 

The  use  of  the  statistical  method  in  solving  a  problem  will  now  be  demonstrated  by  predicting 
the  diffusion  of  acetone  vapor  through  a  soil  profile. 

As  shown  by  the  experimental  and  theoretical  work  described  earlier,  the  migration  of  volatile 
chemicals  in  the  gas  phase  of  unsaturated  soil  depends  primarily  on  soil  water  content,  void 
porosity,  and  the  degree  of  interaction  of  the  volatile  chemical  with  soil  water  and  minerals.  If 
the  interaction  processes  proceed  very  rapidly  compared  with  the  diffusion  process,  local  equilib¬ 
rium  can  be  assumed  to  exist  between  the  components  in  the  gas  phase,  and  liquid  and/or  solid 
phases.  In  the  case  where  a  linear  relationship  for  partition  of  component  between  the  vapor  and 
immobile  phases  holds,  the  migration  of  volatile  chemicals  can  be  described*  by  eq  92. 

In  determining  the  adsorption  of  acetone  by  Laguna  Joyuda  soil,  it  was  found  that  the  following 
simple  relation  holds  when  the  concentration  of  acetone  is  low: 


S  kCj  (18) 

where 

C  --  concentration  of  acetone  per  unit  volume  of  the  soil  gas  phase  (g  cm'J) 

S  concentration  of  acetone  per  unit  volume  of  interfacial  water  and  minerals  (g  cm's) 

k  =  adsorption  coefficient  of  acetone  by  earth  media. 

If  physical  variables  such  as  the  water  content  do  not  change  significantly  during  the  time 
period  under  consideration,  the  migration  of  acetone  can  be  described  mathematically  in  Cartesian 
coordinates  as 


dCx 
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where 
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The  value  of  DR  can  be  estimated  using  the 
relationship 

Dh  0.66*  D2.  (73) 

Combining  eq  109  and  73  one  obtains 

0.66  Do 

DAx)  - 1  ,(x).  (119) 

k(x)  »  1 

The  migration  of  acetone  in  soil  from  the 
Laguna  detection  site  for  the  situation  shown  in 
Figure  10  will  now  be  analyzed.  In  this  case,  the 
soil  profile  consists  of  two  distinct  layers,  namely 
a  moist  and  less  compact  layer  overlying  a  dry, 
compact  horizon.  It  is  assumed  that  an  instanta¬ 
neous  line  source  of  acetone  is  deposited  at  a 
depth  of  20  cm  from  the  surface.  The  computation 
on  the  migration  of  acetone  is  performed  as  follows. 

For  the  sake  of  convenience,  the  following 
dimensionless  variables  are  introduced: 


Figure  10.  Model  of  soil  profile  indicating 
position  of  line  source  of  acetone. 
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(113) 


(114) 


Since  D  depends  only  on  the  coordinate  (Xj)  eq  108  reduces  to 


ft  _  iA 

dt+  2^f  dxf 


[2Dj(Xj)Cj]  -  _ 

dx\ 


dD+(x|)^ 

dx\  ) 


(115) 


The  problem  can  be  described  as  a  corresponding  diffusion  process  <f  t(£ lt.  £2t)  with  a  diffusion 
matrix 


B 


(116) 
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and  a  transition  vector  a  l(dD|  dx|),  0,  0|.  From  eq  98  one  obtains  a  matrix  R  as 


Finally,  the  following  formula  is  obtained  for  flt  and 


dD\ 


t  — -  )  \  t  f  2D\  Aij.(m) 

It  mi  1  it  m  \  ^  f  /  t  m  1 


i.  + 

2t f  m 


2D[  Aij2(m) 


tu  —  0,  1 , , ,  i  n  1, 


(117) 


(118) 

(119) 


Since  Dj(x|)  has  discontinuity  at  the  interface  of  the  two  soil  layers,  it  is  obvious  that  the 
Lipschitz  condition  is  not  satisfied.  However,  it  is  possible  to  approximate  the  function  D  j(x+) 
as  close  as  one  wishes  by  another  function  D+(x+)  which  satisfies  all  requirements  without 
distorting  the  physical  picture.  Generally,  a  sufficient  condition  for  a  function  f(x+)  to  satisfy  the 
Lipschitz  condition  is  that  f(x+)  has  bounded  first  derivatives.  In  this  case,  it  is  easy  to  see  that 
if  Drf(x*)  has  bounded  second  derivatives,  eq  118  and  119  have  a  unique  solution.  A  discontinuity 
at  x  jf  0  was  connected  smoothly  by  a  polynomial  of  the  third  order  as  follows: 
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xQ  <  x|  <  1.5. 


(120) 


It  was  assumed  that  the  concentration  of  acetone  above  the  soil /atmosphere  interface  is  virtually 
zero.  In  the  theory  of  probability  such  a  boundary  condition  is  referred  to  as  an  absorbing  barrier. 

If  the  path  of  £t  goes  beyond  the  boundary,  there  is  no  chance  for  the  path  to  reenter  the  boundary. 

The  calculated  distribution  of  acetone  in  the  soil  after  two  days  is  indicated  in  Figure  11  where 
the  final  points  of  each  sample  path  are  plotted.  The  point  density  is  proportional  to  the  acetone 
concentration.  In  this  computation,  the  time  increment  was  taken  as  15  minutes  and  x0  as  0.1.  The 
proper  choice  of  xQ  depends  on  the  value  of  the  time  increment.  xQ  should  be  selected  so  that  the 
sample  path  will  include  at  least  one  point  in  theregions  (-xQ  <  Xj  <  x0)  whenever  it  traverses 
the  boundary,  xj  =  0.  Now,  in  an  idealized  situation  in  which  acetone  is  deposited  as  an 
instantaneous  line  source  in  an  infinite,  homogeneous  soil  having  diffusion  coefficient  Dt,  the 
concentration  distribution  would  be  given  by 
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Figure  11.  Distribution  of  acetone  after  2  days. 

from  which  it  is  easy  to  see  that  the  distribution  of  acetone  would  be  symmetrical.  This  is  also  the 
situation  which  exists  in  the  present  case  before  the  concentration  distribution  is  distorted  due  to 
the  effect  of  soil  inhomogeneity  and  the  atmosphere/soil  interface.  Although  the  distribution  is 
symmetrical  with  respect  to  the  plane  x2  0,  eventually  the  lower  soil  horizon  acts  as  a  barrier 
(low  diffusion  coefficient)  retarding  diffusion  from  the  upper  layer.  Thus,  acetone  accumulates  in 
the  vicinity  of  the  interface.  The  opposite  situation  exists  at  the  atmosphere/ soil  interface  due  to 
the  relatively  rapid  diffusivity  of  acetone  in  air. 

The  method  based  on  stochastic  diffusion  processes  can  be  applied  successfully  to  much  more 
complicated  problems  than  the  example  discussed.  For  instance  the  properties  of  soil  change  not 
only  with  location  but  also  with  time  due  to  the  influence  of  climatic  variation.  Once  chemical 
effluvia  unique  to  buried  explosives  are  characterized,  and  if  proper  geological  and  meteorological 
data  are  both  available  it  should  be  possible  to  use  this  method  to  predict  the  appearance  of  effluvia 
in  the  atmosphere  under  different  natural  environmental  situations. 
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